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Abstract 

A lattice calculation of the form factors that determine the "hadron- 
ization ratios", such as Rk* and R^f,, where Rk* = [^{B K*^)/T{b — 
S'j)], is presented in the quenched approximation. Lattice data shows 
strong evidence for the scaling law suggested by heavy quark symme- 
try for one of the form factors {i.e. T2). The data also gives strong 
support for the simple pole ansatz for the dependence of T2 in the 
range of available q^. We thus find T2(0) = .10 ± .01 ± .03 yielding 
Rk- = (6.0 ± 1.2 ± 3.4)%; we also find R^ = (6.6 ± 1.3 ± 3.7)%. 



The loop decays of the 6-quark have long been noted for their capacity 
to provide important tests of the Standard Model (SM) Since many of 
these decays are short distance dominated at the quark level, their inclusive 
rates are amenable to perturbation theory. Thus for inclusive processes re- 
liable predictions can be made. The fact that the 6-quark has a relatively 
long lifetime facilitates experimental tests of the theoretical predictions. In 
particular, in the SM, the simple decay 6 — > s + 7 is predicted to have a 
branching ratio which varies from ^ 2 x 10^^ to ^ 4 x 10~^ as the top quark 
mass varies from 100 to 200 GeV 0. Assuming three generation unitarity, 
along with Vcs — Vtb — I, it is easy to see that 6 ^ s + 7 is independent of 
CKM angles to a very good approximation. Furthermore, 6 ^ s + 7 is 
also noted for its sensitivity to extensions of the SM . Finally, we note that 
b — » s{d) + 7 can also lead to a determination of Vts and Vtd- However, the 
full potential of the loop decays of the 6-quark is very difficult to capitalize 
upon unless reliable theoretical predictions can also be made at the exclusive 
level. 

Indeed, the inclusive process b s + j is challenging to measure exper- 
imentally; whereas a corresponding exclusive mode [i.e., B ^ K* +7) has 
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a distinctive signature and is much more accessible to experiment. Thus 
a meaningful confrontation between experiment and the underlying elec- 
troweak theory can be facilitated through a knowledge of the "hadronization 
ratio", Rk*'- 

which is the probability for the formation of the K* . The evaluation of this 
ratio by continuum methods has proven to be extremely difficult. This is 
reflected in the wide range ~ 1 to ~ 97% in the value of Rk*, as calculated 
by quark models, QCD sum rules, heavy quark symmetry (HQS) extended to 
include the s-quark, etc. (see Table 1) 0. Under the circumstances, the hard 
earned experimental determination of the branching ratio for B —>■ K* + 7 
may only be used to select amongst various models of hadronization rather 
than to test the underlying theory. It is thus clearly important to explore 
the use of lattice methods for treating such exclusive decays of B mesons. 
At the quark level the decay is described by an effective Hamiltonian 
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Heff = Gg^g^g^{mu^) Vfs s ( X ) (T^^ ( X ) F^'^ix) (2) 

where F^'^ is the photon field strength tensor, ha = ^-^^b, and the c number 
coefficient Gg^g^g.j^ijnt, n) depends on all the three gauge couplings of the 
Standard Model (SM), the mass (mj) of the top quark and a renormalization 
point 11. The dependence on the CKM angle Vts has also been factored 
out. Vtb is assumed to be 1, and the negligibly small u quark contribution is 
ignored. 

As usual [|] the lattice is used for a non-perturbative evaluation of the 
matrix element = (y{k)\J^\P{p))^ where P is the initial pseudoscalar 
heavy-light meson, V is the final vector meson, = sa^^q^jhn = (f + a)^ is 
the current, with and the vector and axial parts, and q=p — k\s the 
4-momentum of the photon. In general, the Euclidean matrix element can 
be parameterized in terms of three form factors : 



-7] ■ q{p + k)^]T2{q^) + r] ■ q 



rrifj — m 



iP + k), 



Tsiq') (3) 



V 
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(Our 7 matrices obey {jij.,ju} = "^Sf^u, and momenta are defined by — 
{E,ip), so that — on shell.) For a lattice calculation it is simpler to 
note that the term arises purely from the vector piece of and the 

T2 and T3 terms given above arise from the axial piece. The third term does 
not contribute when the photon is on shell. Furthermore, at the end-point, 
(g^ — q^ax = {''^H — where the final and initial mesons are both at 

rest, T3 term does not contribute to the axial matrix element. Since also 
at that kinematic point no momentum injection is required, 72 (g^^^,) can be 
readily, and rather cleanly, evaluated on the lattice. Although g^ = (the 
point of direct physical interest) is not exactly accessible to the lattice, in 
many instances the parameters used in the current simulation do allow to 
be extremely small, i.e. /rri^j < .1. Finally we note that using the identity 
cr^y 75 = — |e^i,ApCr'^^ one can show that 

T2(0) = Ti(0). (4) 

Now the hadronization ratio of interest takes the simple form (for rus <S 

rrib): 
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Ti(0)p . 



(5) 



On current lattices = (or near that point) is inaccessible for very heavy 
meson masses, say tjih > 3.5 GeV. So at ttih ~ ms, Ti{0) is not di- 
rectly calculable. However HQS allows one to predict the behavior of 
'^^{(Imax) large uih- Indeed, when = g^^^ no large momenta is trans- 
ferred to the recoiling light hadron, so a straightforward argument shows 
that y/^HT2{qmax) ~^ coHst. (up to logarithms) as rriH oo. This makes 
possible a controlled extrapolation of 22(5^^3.). strategy on the lattice 
will thus be: 

1) test pole dominance of T2 at fixed rriH, to the extent that the data 
allow, by deducing T2(0) from T2{ql^^^) using the equation: 



and comparing to Ti(0) using eq. (4). Ti(0) is also obtained using pole 
dominance, but only from Ti at small values of {(f/^^ < 0.1). Pole 
dominance does not appear to work well for Ti(g^) with large g^. 

2) extract 72(g^^^) at ttih = by fitting the data to the form suggested 



T2(0)=T2(gLJ 1 




(6) 
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by HQS, namely: 



( 



1 



) 



(7) 



rriH 



3) use pole dominance for T2 at itih = ms to deduce Ti(0) = T2(0) from 



We remark that in testing pole dominance, we have simply used the 
pseudoscalar mass in eq.(6) because in the limit of large mn, HQS implies 



that step 3) uses pole dominance over a wider range of {imaxl^'^'B ~ 0.65) 
than can be explicitly checked in step 1) {qmax/^'H ~ 0.3). We attempt to 
estimate below the systematic error associated with this step. 

We mention the following technical points, regarding the lattice calcula- 
tions, in brief 

1) The recently proposed normalization of the Wilson quarks on the lattice 




that resonances of different spin parities become degenerate p[. Note also 



M MM- 




(8) 



where 
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am = In 



1 

2k " . 



and K = k,/8kc {hc is the critical hopping parameter) is used. Thus the 
leading corrections that become important as am gets large are automatically 
included. 

2) For the renormalization of the tensor current we incorporate the correction 



calculated in lattice weak coupling perturbation theory to one loop order [12 . 



However, following Lepage and MacKenzie |]rT|, the tadpole contribution 
is removed from the correction (it is already included in eq. (8)), and a 
"boosted" value of = gv{l/a) is used. 

We have done the calculation of Ti and T2 on four different sets of lattices: 
A) p = 6.3, X 61 (20 configurations, a~^ = 3.01 GeV), Bl) P = 6.0, 
243 X 39 (8 configurations, a~^ = 2.29 GeV), B2) p = 6.0, 24^ x 39 (a 
second set of 8 configurations, a~^ = 2.29 GeV), and C) /3 = 6.0, 16^ x 39 
(19 configurations, a^^ = 2.10 GeV). The lattice spacings given above are 
determined through a calculation of /^r with the same point sources that are 



used here |12]. The "i?" is always taken at rest; the "i^*" is given lattice 



momentum (0,0,0), (1,0,0), (1,1,0) or (2,0,0), with (2,0,0) used only on 
Bl and B2. Preliminary results of this computation have been presented 
previously. ||14 



We first work in the case when the masses of the two light quarks are held 
equal. Experimentally this situation corresponds to the decay, for example, 
i?s ^ </) + 7 . For the light quark we use k = .148 at /? = 6.3 and k = .152 
aX (3 = 6.0, yielding a vector meson in the final state with mass ~ 1.3 
GeV. The dependence of the amplitude on the heavy quark mass is then 
studied. Specifically, for (3 = 6.3, we use k = 0.140, 0.125, 0.110 and 0.100 
for the heavy quark. Results are given in Table 2; the last column shows 
that ^/rrlH T2{q^ax) approximately constant. We then fit the data to the 
two parameter form (eq. (6)) suggested by HQS taking the correlations in 
the data into account through covariant fits. For the (3 = 6.3 data we find 

= .806 ± .069 (GeV)5, A2 = -.545 ± .082 (GeV)i, {x^dof ^ 2.3/2). 
Thus 

T2{mH = ruB^^q'^ax) = -304 ± .030 . 

We now discuss the systematic errors on T2, first considering those rele- 
vant to i?s ^ + 7. To correct for the physical s-quark we also study the 



matrix elements with k = 0.149 (corresponding to vector meson of mass 
about 1.1 GeV), aX (3 — 6.3. We find a shift in T2, from its value k — 0.148 
of about —7%. Extrapolating to the physical s-quark would give a shift of 
— 10%. In passing we mention that a similar study of our lattices at j3 = 6.0 
indicates a smaller error than the 10% seen at /3 = 6.3. 

We now assess the systematic error due to the use of heavy quarks with 
am > 1. For that purpose, we fit to the two parameter form using the two 
lightest quarks from our heavy set (of four) at (3 — 6.3; i.e. we retain only 
K = 0.140 and k = 0.125. We find a shift in the value of T2 of 3.1%. 

To estimate scale breaking errors we compare the fit for the (3 = 6.3 data 
with the heavy quarks a,t k, — .140 and .125 to the fit for the P — 6.0 data 
with the corresponding heavy quarks at k — .135 and .118. We attribute the 
difference of 12.2% to scaling violations. 

The systematic errors due to finite size effects are deduced by comparing 
the value of T2 on our 16^ lattice with the one on the 24^ lattice, both at 
p = 6.0. We find a difference of 9.4%. 

Adding in quadrature the errors due to the four sources mentioned above 
we find a total systematic error of 19%. In passing we note, however, that 
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the systematic error due to each of these four sources is actually smaller than 
the statistical errors in the appropriate subset of data. It is, therefore, quite 
hkely that the estimate of 19% is a conservative one. Thus, we arrive at: 



T2{mH = mB^,qLx) = -304 ± .030 ± .057 . 



(9) 



Table 3 summarizes our test of the pole dominance for T2. By examining 
the agreement between 72(0) and Ti(0) we see that, within the available range 
of Qmax/'^'H — 0-3, the pole-model seems to work very well. We must note, 
however, that in the actual physical reactions of interest (i^axl^% approaches 
about 0.65. To estimate the error involved, we note that the biggest difference 
between Ti(0) and T2(0) is ~ 9% (for lattice C). Scaling by the increased 
range in (f for the physical reaction (0.3 — > 0.65), we arrive at an error of 
20%. Since the data points with higher q^ax/'^'n Table 3 seem to support 
pole ansatz just as well as those with lower values of qmaxl^^Hi this is likely 
to be an overestimate, but we wish to be conservative. Using eqs.(4), (5), 
(6) and (9) we thus find: 



^^^(0) = 7f ^^^(0) = 0.104 ± 0.010 ± 0.028 . 



(10) 
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= (6.6± 1.3±3.7)% , (11) 

which is the hadronization ratio for i?^ — > + 7- Note that in this calculation 
we have taken rrih = 4.5 GeV, so that we may use the result for BR(6 — >■ 
sj) given by Misiak 0. A 13% uncertainty is added in quadrature to the 
systematic errors on R corresponding to an assumed 200 MeV uncertainty 
in mf,. 

Next we turn our attention to B ^ K* + 7. For this purpose we study 
matrix elements with unequal masses for the light quarks. For example, at 
/3 = 6.0 we take the spectator quark with k = 0.154 and the "s" quark 
(corresponding to the light quark that results from the weak decay of the 
6-quark) with k = 0.152. Furthermore, we have to extrapolate in the masses 
of these two quarks. In particular, the spectator quark requires extrapolation 
to the chiral limit {i.e. Kc = 0.157 at /? = 6.0). For this study we use the 
P = 6.0, 24^ lattice as it has the largest physical volume. This lattice has two 
independent sets of configurations with eight configurations in each sample. 
As a result we find that extrapolation to the physical limit causes a shift 
in T2 from its value calculated with degenerate light quarks ( k = 0.152 at 
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P = 6.0) of about +7%. We shift the central value of T2 accordingly and also 
include this additional 7% in the systematic errors. Consequently we arrive 
at: 

T2{mH = TUB, ql,ax) = -325 ± .033 ± .065 . (12) 
Once again we use pole dominance to get: 

T^^^'{0) = T^^^\Q) = 0.101 ± 0.010 ± 0.028 , (13) 

Rk' = {6.0 ±1.2 ±3 A)% . (14) 

Now, as mentioned earlier, the inclusive branching ratio for 6 — > S7 is 
predicted to he in the range of about (2 — 4) x 10~^ depending on m^. Thus, 
in the SM, there is a bound, BR{b — > s + 7) < 4 x 10"^, corresponding to 
rrit ~ 200 GeV. Combining this upper bound with the above lattice result 
one gets: 

BR{B K*-f) < (2.4 ± 0.5 ± 1.4) x 10"^ . (15) 
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We recall now the recent CLEO result [|I5|1: 



BR{B K*-i) = (4.5 ± 1.5 ± 0.9) x 10"^ . (16) 

Given the size of the errors in the lattice calculation, as well as in the 
experiment, the CLEO result is certainly not inconsistent with the expec- 
tations based on the lattice. However, we note that the numbers seem to 
mildly favor a rather heavy top quark. 

In an attempt to quantify the statement about mt we note that the ex- 
perimental result along with the lattice implies: 

BR{b S7) = (4.5±1.5±0.9)xl0"V(6-0±1.2±3.4)xl0~2 ^ (7.5±5.4)xl0~^ 

where we have assigned a ~ 70% combined error to the lattice plus the 
experimental result. At the l-cr level one then finds nit ^ lOOGeV^. However 
modest improvements in the the lattice and/or experimental results could 
produce a rather stringent bound. 

To summarize, we have used lattice methods to evaluate the form factors 
for the radiative B transitions. We emphasize that the heavy quark limit 
of QCD P) enters in an important way in making this calculation feasible 
on current lattices. We also want to highlight two drawbacks of the present 
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effort. First, numerical limitations did not allow us to check pole dominance 
for the specific value of the momentum transfer relevant to the experiment. 
However, lattice parameters did allow us to check the pole model for T2(0) for 
an appreciable range of momenta, giving strong support to its validity. We 
have included what we believe is a conservative estimate of 20% systematic 
error due to the use of pole dominance. The second limitation is, of course, 
the quenched approximation. It is generally believed that with the use of a 



physical quantity {e.g. /,r in our work ||T2[) to set the scale for the lattice 
calculations, errors due to quenching are hkely to be quite small, perhaps 
< 10%, in the form factors of interest here. It is, therefore, unlikely that the 
present limitations would seriously affect our results, given the ~ 28% error 
in amplitude. Quenched simulations are now in progress that should allow 
us to improve the calculations to the 10-15% level. At that stage errors due 
to quenching may also start to become relevant. 
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Table 1: A sample compilation of the predictions for Rk* = 



[r{B - K*^)/T{b s-f)]. See Ref. 5. 

Author (s) Rk* 

O'Donnell (1986) 97% 

Deshpande et al. (1988) 6% 

Domingues et al (1988) 28 ± 11% 

Altomari (1988) 4.5% 

Deshpande et al. (1989) 6-14% 

Ahev et al. (1990) 39% 

Ah et al. (1991) 28-40% 

Du et al. (1992) 69% 

Faustov et al. (1992) 6.5% 

El-Hassan et al. (1992) - 0.7% - 12% 

O'DonneU et al. (1993) ~ 10% 

Ah et al. (1993) 13 ± 3% 

BaU (1994) 20 ± 6% 

This work (6.0 ± 1.2 ± 3. 
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Table 2: Lattice data on four sets of lattices, ki represents 
the intial heavy quark undergoing weak decay, K2 the light quark 
emerging from the weak decay. The spectator quark is taken to 
have K2 as well, mn and my are the masses of the initial and the 
final 0~ and 1"" mesons respectively and Vmax = femaa:/"^//]- 



(3{a-^/GeV) 
{'.} 


Lattice 

Sot 


Kl, K2 


airiH 


amy 




T2 {(Irruu') 


^mHl2[.qmax) 

(GeV)^ 


6.3(3.01) 


A 


140,148 


.590 


.422 


.081 


.406 ± 


.046 


.54 ± .06 


{0.151} 




















A 


125,148 


.934 


.422 


.301 


.384 ± 


.044 


.64 ± .08 




A 


110,148 


1.248 


.422 


.443 


.364 ± 


.048 


.71± .10 




A 


100,148 


1.465 


.422 


.508 


.346 ± 


.055 


.73± .12 


6.0(2.29) 


Bl 


135,152 


.894 


.561 


.139 


.409 ± 


.090 


.58± .13 


{0.157} 




















Bl 


118,152 


1.244 


.561 


.301 


.371 ± 


.105 


.63± .18 


6.0(2.29) 


B2 


135,152 


.891 


.566 


.139 


.478 ± 


.090 


.69± .13 




B2 


118,152 


1.245 


.566 


.301 


.415 ± 


.105 


.71± .17 


6.0(2.10) 


C 


142,152 


.734 


.564 


.053 


.470 ± 


.062 


.58± .08 




c 


135,152 


.888 


.564 


.133 


.459 ± 


.065 


.63 ± .09 




c 


118,152 


1.241 


.564 


.298 


.414 ± 


.089 


.67± .14 
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Table 3: Test of the pole model for the dependence of the 
form factors; in particular, that of 72- ^2(0) and Ti(0) are deduced, 
from Ti{q^) and T2{q^) seen on the lattice, by using pole dominance 
i.e. eq. (6). Note r = q^ /m\j- 



Lattice 
Set 




r 


'^max 


TM 


') 




TM 


^2(0) 


A 


125,148 


.002 


.300 


.259 ± 


.035 


.384 ± 


.044 


.259 ± 


.035 


.269 ± 


.032 


Bl 


135,152 


.009 


.139 


.391 ± 


.069 


.409 ± 


.090 


.388 ± 


.068 


.352 ± 


.077 


B2 


135,152 


.009 


.139 


.436 ± 


.092 


.478 ± 


.090 


.432 d= 


.091 


.411 ± 


.078 


Bl 


118,152 


-.034 


.301 


.264 ± 


.050 


.371 ± 


.105 


.272 ± 


.051 


.260 ± 


.073 


B2 


118,152 


-.034 


.301 


.316 ± 


.110 


.415 ± 


.100 


.327 ± 


.113 


.290 ± 


.070 


C 


118,152 


-.069 


.298 


.300 ± 


.039 


.414 ± 


.089 


.321 ± 


.042 


.291 ± 


.062 
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